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ABSTRACT 

The nature of solar wind (SW) turbulence below the proton gyroscale is a topic that is being investigated extensively 
nowadays, both theoretically and observationally. Although recent observations gave evidence of the dominance 
of kinetic Alfven waves (KAWs) at sub-ion scales with a> < co C i, other studies suggest that the KAW mode cannot 
carry the turbulence cascade down to electron scales and that the whistler mode (i.e., co > &> ci ) is more relevant. 

Here, we study key properties of the short-wavelength plasma modes under limited, but realistic, SW conditions, 
typically /3, > fi e ~ 1 and for high oblique angles of propagation 80° $7 0kB < 90° as observed from the Cluster 
spacecraft data. The linear properties of the plasma modes under these conditions are poorly known, which contrasts 
with the well-documented cold plasma limit and/or moderate oblique angles of propagation (0kB < 80°). Based on 
linear solutions of the Vlasov kinetic theory, we discuss the relevance of each plasma mode (fast, Bernstein, KAW, 
whistler) in carrying the energy cascade down to electron scales. We show, in particular, that the shear Alfven mode 
(known in the magnetohydrodynamic limit) extends at scales kp, > 1 to frequencies either larger or smaller than 
co c i, depending on the anisotropy k\\/k±. This extension into small scales is more readily called whistler (co > co c \ ) 
or KAW (co < &> c i), although the mode is essentially the same. This contrasts with the well-accepted idea that the 
whistler branch always develops as a continuation at high frequencies of the fast magnetosonic mode. We show, 
furthermore, that the whistler branch is more damped than the KAW one, which makes the latter the more relevant 
candidate to carry the energy cascade down to electron scales. We discuss how these new findings may facilitate 
resolution of the controversy concerning the nature of the small-scale turbulence, and we discuss the implications 
for present and future spacecraft wave measurements in the SW. 

Key words: magnetohydrodynamics (MHD) - plasmas - solar wind - turbulence - waves 
Online-only material: color figures 


1. INTRODUCTION 

The nature of the solar wind (SW) turbulence below the ion 
scale (typically p ,■ ~ 100 km, corresponding to an observed 
frequency in the spacecraft frame of / sc ~ 0.5 Hz) has attracted 
considerable interest in the space and astrophysical communities 
in recent years. This has been encouraged in particular by 
recent observations from the Cluster mission that provided 
the most complete and detailed picture of the SW turbulence 
cascade from magnetohydrodynamic (MHD) scales ( L P> p, ) 
to electron scales (L ~ p e ~ 1 km; Sahraoui et al. 2009, 
2010b; Kiyani et al. 2009; Alexandrova et al. 2009; Chen et al. 
2010). Determining the nature and properties (e.g., scaling, 
anisotropy) of the turbulence at small scales is indeed a crucial 
point to understanding the problems of energy dissipation and 
heating, particle acceleration, and magnetic reconnection in 
space and astrophysical plasmas (Schekochihin et al. 2009). 
Recent Cluster observations provided clear evidence that SW 
turbulence cascades below the ion scale p ,■ down to the electron 
scale p e where dissipation becomes important and the spectra 
steepen to ~/ sc _ “, with a > 4 (Sahraoui et al. 2009, 2010b). 
The spectrum thus formed was termed the dissipation range 
while the range of scales between p, and p e has been termed 
the dispersive range in reference to the dispersive nature of the 
plasma modes at those scales (Stawicki et al. 2001). In a case 
study, Sahraoui et al. (2010b) performed a detailed analysis 
of the energy cascade from MHD to sub-ion scales using a 
multipoint measurement technique called A:- filtering (Pin 5 on & 
Lefeuvre 1991; Sahraoui et al. 2003b, 2010a; Naritaet al. 2010a, 
2010b; Tjulin et al. 2005). The results clearly showed that; 


(1) the magnetic turbulence is strongly anisotropic (k± k y) 
down to the observed scale k±p, ~ 2; (2) the cascade is 
consistent with kinetic Alfven wave (KAW) turbulence as 
proposed in Howes (2008) and Schekochihin et al. (2009) with 
frequencies in the plasma rest frame co < 0. 1 oi CI , although the 
frequencies in the spacecraft frame reached 20tt> c ;; and (3) the 
turbulence undergoes a transition range near the ion scale p, 
characterized by a steepening of the spectrum from V/ 1 ' 6 to 
k . 4 ' 5 , which has been interpreted as due to Landau damping 
of magnetic energy into ion heating (Howes et al. 2008). The 
remaining energy cascades following a power law ~/“ 28 
down to the electron scale where the energy was suggested 
to dissipate into electron heating. The wavenumber spectra at 
those frequencies could not have been directly measured using 
the ^-filtering technique due to the limitation imposed by the 
satellites separations, which are larger than 100 km (Sahraoui 
et al. 2010c). However, one-dimensional wavenumber spectra 
can be inferred indirectly by means of the Taylor frozen-in- 
flow approximation provided that it is valid (Bale et al. 2005; 
Sahraoui et al. 2009; Alexandrova et al. 2009). Similar analyses 
of magnetic turbulence at MHD and sub-ion scales using the 
^-filtering technique have been carried out recently and have 
confirmed most of these results (Narita et al. 2010a, 2010b, 
201 1 ; F. Sahraoui et al. 2011, in preparation). 3 


3 Narita et al. (201 1) have claimed to have observed several plasma modes, 

other than the KAW, based on the random spread of the observed dispersion 

relations. However, given the absence of any statistical error bars in those 

observations, one cannot reliably interpret those results. A detailed discussion 
of these important issues will be given elsewhere. 
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From the theoretical point of view, the debate is generally 
polarized between those who advocate for the relevance of KAW 
turbulence to explain observations of small-scale SW turbulence 
(Howes 2008; Howes et al. 2008, 2011b; Schekochihin et al. 
2009) and those who believe in the necessity that whistler (or 
another type of) turbulence should exist at small scales (Stawicki 
et al. 2001; Gary & Smith 2009; Podesta et al. 2010). Indeed, 
Podesta et al. (2010) have used a toy model of energy cascade 
where the linear damping rate of the KAW waves is considered 
and showed that the KAW cascade cannot reach the electron 
scale because the energy flux of the cascade vanishes at scales 
kpi ~ 20. However, recent Gyro Kinetic simulations, which self- 
consistently contain kinetic damping of the plasma modes, did 
not confirm that conclusion and showed rather that the KAW 
mode can carry turbulence cascade down to electron scales 
(Howes et al. 201 lb). Recent two-dimensional PIC simulations 
also showed the formation of a power-law spectrum L[ 5 ' 8 below 
the electron inertial length (Camporeale & Burgess 2011). In 
the framework of incompressible electron-MHD turbulence, 
Mayrand & Galtier (2010) showed that the magnetic energy 
spectra should follow a power law k~ n ^ at scales smaller than 
d e (the electron inertial length), which they proposed to explain 
the steepening of the spectra near the electron scale reported 
in Sahraoui et al. (2009, 2010b). Note, however, that this fluid 
model is non-dissipative and does not consider any damping of 
the turbulence via kinetic effects, which are important in the 
dispersive and the dissipation ranges. 

Here we derive and analyze the linear solutions of the 
Vlasov-Maxwell equations to shed light on new aspects of 
small-scale kinetic plasma modes. Based on the obtained damp- 
ing rates of the waves and their magnetic compressibilities, we 
discuss the relevance of each mode to carry the energy cascade 
in the SW down to electron scales. Unlike the previous papers 
that had dealt with this problem, in this article we focus only on 
the properties of the plasma modes under a restricted range of 
SW parameters, namely, hot plasma with /!,■ > /3 e ~ 1 and high 
oblique angles of propagation, typically 80° ^ 0 rb < 90°, as 
observed recently using the ^-filtering technique (Sahraoui et al. 
2010b; Narita et al. 2011; F. Sahraoui et al. 2011, in prepara- 
tion). This range of parameters has not been addressed carefully 
so far. Filling this gap should help resolve part of the ongoing 
controversies on the problem of energy cascade and dissipation 
in the dispersive range. 

Before studying in detail the plasma kinetic modes, in the 
next section we will first derive the linear solutions of the two- 
fluid theory in the cold and hot plasma cases. Although one 
may question the validity of the fluid description of hot plasmas 
(we do not address that issue here; see Krauss-Varban et al. 
1994; Howes 2009), we believe nevertheless that deriving the 
fluid solutions is a necessary step to better understand the more 
complex solutions given by the kinetic theory. This complexity 
is essentially due to the damping of the modes by kinetic effects, 
which makes them hard to trace at all scales, in addition to 
the rising of several mode conversions (Krauss-Varban et al. 
1994; Li & Habbal 2001). As we will show, the reduced-two- 
fluid theory indeed helps track the various plasma modes and 
understand the important changes in their properties depending 
on the plasma f J > and the obliquity of the waves. It also allows 
us to derive analytically various asymptotic properties of the 
waves that will be used to discuss new properties of the Alfven 
and the whistler modes. In the last section we discuss the results 
and their implications on spacecraft measurements of waves in 
the SW. 


2. LINEAR SOLUTIONS OF THE 
REDUCED-TWO-FLUID MODEL 

We use the reduced-two-fluid theory as described in Sahraoui 
et al. (2003a), whose equations are the following: 


+ V • (n s \ s ) = 0 (1) 

V P s q s 

3/V, s + v s ■ Vv 5 = + — (E + v s x B) (2) 

m s ii s m s 

Ps 

d, — =0 (3) 

( m s n s )Y° 

3,B = —V x E (4) 

V x B = /r 0 J, (5) 


where the subscript .y denotes ions and electrons, y s is the 
polytropic index, and J = ]C 5 q s n s \ s is the electric current. 
The system of Equations ( I )— (5) is essentially equivalent to 
the classical two-fluid model when the non-relativistic and 
quasi-neutrality assumptions are considered. This translates into 
neglecting both the longitudinal and the transverse components 
of the displacement current, which yields the simplified form of 
Equation (5). This system has the advantage of ruling out the 
three high-frequency modes of the two-fluid theory: the optic 
modes (with ~ c, the speed of light) and the Langmuir mode 
o) ~ o> pe . However, it still includes the electron inertia, which 
allows us to describe properly the whistler mode. The system has 
therefore three eigenmodes whose dispersion relations can be 
derived analytically. The main simplification with respect to the 
complete kinetic system that will be studied in the next section 
concerns the introduction of polytropic closure equations that 
excludes any resonant phenomenon, in particular any effect of 
wave damping. 

To derive the linear solutions we linearize Equations ( 1 )— (5) 
by assuming that n s = no + 8n s , P s = Pq + 8P s ,Ji = 
Bo + <5B. If, furthermore, we assume that the perturbations 
vary as e -/( <8f - k r ) ; 0 ne can obtain (after some calculations) the 
following set of equations: 


where 


<5E = BqM,- • 8\i 


<SE = BoM f • 8\ e 


1 = 

8\j — 8v e = — M ■ 5E, 
Bo 


A/I . o) - . k 2 Vthi 2 A 

M,- = - j — I + JYi e k e k + G 

COCOr'l 


co = k 2 V th f 2 = 

M„ = j — I - jy e e k e k + G 

COce CO (Ore 
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(7) 
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Figure 1 . Dispersion relations of the reduced-two-fluid model in the hot (red) and cold (black) plasmas with 7/ = T e , yt = 5/3, and y e = 1 at different angles of 
propagation. The vertical dashed lines correspond to the electron inertial length kd e = 1 for the two values of /?, . The inset in the top panel is a zoom on the very small 
scales of the plot to show the asymptotes co+ and co - . 

(A color version of this figure is available in the online journal.) 


where I is the identity tensor, G = e y e x — e x e v , and e k is the 
unit wave vector assumed to be in the XZ plane 

/ sin 0 k B \ 

e k = k/k = 0 . 

\cos 0 kB / 

Bo = Boe z , co cs = eB/m s , Vthj, and V As are the back- 
ground magnetic field, the gyropulsation, the thermal, and 
the Alfven speeds of the particle s, respectively. Regrouping 
Equations (9)-(l 1) we obtain the final equation 

[m,- ' - M e ' - m] • <$E = 0, 

which admits solutions only if 

Det^M, - M, - M e • M • M, j = 0. (12) 


Equation (12) yields the final dispersion equation of the reduced- 
two-fluid model 

c 4 yX 6 -c 2 [(l +n ei )(l+2p iy )X 4 + (l+iu, 2 ei )yX 6 ]Y 2 

+ [A(l + 1 + c 2 + fry)X 2 + [c 2 ( 1 + nl) 

+ l l ei ( 1 + + 2 Ak)]^ 4 + X 6 ]Y 4 ~ (1 + M-ci) 

[A 2 (l + /4) +2/^(1 + !M ei )X 2 + H 2 ei X 4 ]Y 6 = 0, (13) 

where X = kp t , Y = o>/a> ci , p ; = V thi /co ci , fr = V th 2 /V A 2 , 
c = cos0 kB , Pei = nie/m-i, and y = y e + y t is the total 
polytropic index. The solution of Equations (13) are plotted 
in Figures 1 and 2 in low and high plasmas for the given 
angles of propagation and SW parameters. 

At low frequencies (co < &> c j) we found that the slow 
magnetosonic mode in both low and high /3 has an asymptotic 
frequency 

<A) ~ &>ci cos0 kB . 
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Figure 2. Dispersion relations on a log-log scale of the reduced-two-fluid model 
in the hot (red) and cold (black) plasmas showing the connection between high- 
and low-frequency modes (the same description as in Figure 1 applies). 

(A color version of this figure is available in the online journal.) 

At high frequencies {to > oi cl ) both hot and cold plasmas have 
two modes with different asymptotes. Figure 2 shows a similar 
dispersion curve in log-log scale. One can see that the first 
mode is connected at low frequency {to < o> cl ) to the shear 
Alfven mode (known in MHD) and at kd e 'Jp 1 develops an 
asymptote &>_ given by 

co - ~ o> C e cos 0kB- 

It is worth recalling here that the asymptote to_ is usually 
attributed to the whistler mode, which in this case is connected 
at low frequency to the Alfven mode and not to the fast 
magnetosonic mode. Therefore, in the following we will refer 
to this mode as the Alfven-whistler mode for the sake of 
clarity. We recall here that the designation of slow and fast 
modes is based upon the ordering in the phase speeds of the 
modes. While this ordering remains valid in fluid theories at 
low and high frequency, it can be totally violated in kinetic 
theory due to possible interconnections between low- and high- 
frequency modes (Krauss-Varban et al. 1994; Li & Habbal 
2001). Figure 2 also shows that the second high-frequency mode 
is connected at low frequency (both in high and low ft ; ) to 
the fast magnetosonic mode. Therefore, we will refer to this 
mode as the classical-whistler mode. In the low ft; limit the 
classical-whistler mode undergoes significant curvature near 
oo+ at all angles of propagation before converging toward the 
magnetosonic dispersion co ~ kC s , where C s = y,, + y, V,\ u 

is the sound speed and T ei = T e /Tj. From Equation (13) one 
can easily show that in the cold case {ft; = 0) the “asymptote” 
oo+ is given by 

00+ ~ y/ a> ce 2 COS 2 0 kR + CO c iCO ce , 

which has two interesting limits (for co c ; co ce ): 

1. 

lim oo+ ~ lim to- = co ce , 

$kB ^0 #kfi— ^0 

2 . 

lim oo+ = +J co C [CO C e- 

$kB ^ § 

The first limit shows that the branch co+ tends toward 
the known whistler asymptote to+ ~ co ce in quasi-parallel 
propagation as shown in Figure 1 (top panel), while the second 


limit shows that to+ tends toward the lower-hybrid frequency 
oolh ~ y/oOcitOce ~ 42&> c ; for quasi-perpendicular propagation. 
In contrast, and as can be seen in Figure 1 , the asymptote &>_ of 
the Alfven-whistler mode continues decreasing as 0 k B — » tt/ 2 
in both low and high ft; and becomes even smaller than oo c \ at 
all scales for the angles of propagation 

$kB > $crit. — 1 COS (/X e /). 

For the real mass ratio /i,,; = 1 / 1 836 one obtains 0 cnl = 89°97. 
This implies that at quasi-perpendicular propagation (i.e., 0 k B > 
6? C rit. ) there are two asymptotes below ao c p one at extremely low 
frequency, which is the slow magnetosonic mode discussed 
above {ooq ~ oo c \ cos and the other one is the Alfven- 
whistler mode with the asymptote <w_ ~ w ce cos 0 k B • This is 
the first important conclusion of this part. It is worth recalling 
that these two modes were found to be linear solutions of 
the Hall-MHD equations in the incompressible limit (Sahraoui 
et al. 2007; Galtier 2006), and become degenerate at large 
(MHD) scales. They were referred to, respectively, as the 
Alfven and the whistler modes. In such incompressible fluid 
theories, the classical-whistler mode is indeed ruled out by the 
incompressibility assumption (i.e., ft — > oo). All turbulence 
theories built in the framework of those incompressible fluid 
models are thus based on the slow and Alfven-whistler modes 
discussed here (Galtier 2006, 2008). The second important 
finding from Figure 1 is the major change that the classical- 
whistler mode undergoes in high ft ; : when approaching quasi- 
perpendicular propagation the curvature near co+ disappears 
and the dispersion curve tends toward the magnetosonic mode 
co ~ kC s . As we will show in the next section, when kinetic 
effects are considered, the classical-whistler mode does not 
simply extend above oo c [, but rather splits there into different 
ion Bernstein modes. Because the slow magnetosonic mode is 
strongly damped in high ft plasma by kinetic effects, the only 
relevant mode to carry the energy cascade down to the electron 
scale is the new Alfven-whistler mode. 

3. LINEAR SOLUTIONS OF THE VLASOV-MAXWELL 
EQUATIONS FOR HIGH OBLIQUE PROPAGATION 

In this section we solve numerically the linear 
Maxwell- Vlasov equations using the WHAMP code (Ronnmark 
1982), and compare to the previous linear solutions of the re- 
duced hot two-fluid theory. We assume Maxwellian distribution 
functions of electrons and ions and use realistic SW param- 
eters reported in Sahraoui et al. (2009). Part of the obtained 
dispersion relations and damping rates are given in Figure 3 
(the slow magnetosonic mode was found to be heavily damped 
and thus not plotted here, nor is that mode studied in the rest of 
this paper). The plot shows that the fast magnetosonic modes 
at high oblique angles undergo resonances at the ion cyclotron 
frequency, which contrasts with the better known fast mode in 
quasi-parallel propagation or low ft;. This result is confirmed 
by the analysis of the wave polarization. In Figure 4, we plot 
the phase of the electric field component E y whose sign gives 
the sense of polarization: Arg {E y ) > 0 (resp. < 0) for right 
(resp. left) hand polarized waves. We see a clear transition from 
right- to left-hand polarization of the fast mode near kp; ~ 1 
(corresponding to co ~ o) CI ). 

In Figure 3, we also find the kinetic counterpart of the 
Alfven-whistler mode derived from the hot two-fluid theory. The 
Alfven-whistler mode extends the shear Alfven mode to small 
scales as can be seen in the inset of Figure 3 (this result was also 
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kp. 


Figure 3. Linear solutions of the Maxwell-Vlasov equations: dispersion 
relations (blue) and damping rates (red) for the angles of propagation 80° 

$kB ^ 89° . The inset is a log-log plot of the same dispersion relations to show 
the connection between low- and high-frequency modes. 

(A color version of this figure is available in the online journal.) 

found in Figure 1 of Krauss-Varban et al. 1994, although not 
discussed in detail in that paper). The mode becomes dispersive 
at scales kp, > 1 and develops frequencies larger (resp. smaller) 
than (o c ; for t\is < 88° (resp. ^ 88°) up to the scale kpj ^ 10. 
We refer to the branches u> < a> a and co > &) cl , respectively, as 
the KAW and the whistler modes. Note that the limit w = (» c , 
is reached at different spatial scales depending on the value 
of the angle 0 rb- Figure 4 shows that in contrast to the fast 
magnetosonic mode, the Alfven- whistler mode has a right-hand 
polarization at all scales and does not undergo significant change 
near kpt ~ 1. Nevertheless, the Alfven-whistler modes show 
features related to wave-particle resonances near the harmonics 
of ions (because they are not strictly circularly polarized). This 
can be seen in Figure 5 which shows the enhancement of 
the damping of the Alfven-whistler modes near the harmonics 
a) N = Nco c i — k\\ Vth, (a similar observation can be made in 
Figure 9 at other angles of propagation). 



Figure 5. Exponential of the normalized damping rate of the Alfven-whistler at 
# kB = 87° showing enhancement of the damping near the harmonics of ions. 

Figure 6 shows the Alfven-whistler solutions extended to high 
frequencies and small scales. We observe that the damping of the 
Alfven-whistler mode becomes more important when departing 
from ()y\\ ~ 90° toward less oblique angles. For 0 rb = 89.99 
the solution extends down to the electron gyroscale p e , where 
the damping rate remains small, 4 y / co r ~ 0.4. For less oblique 
angles, the Alfven-whistler mode develops frequencies higher 
than op, but they are subject to stronger damping. This can be 
seen clearly in Figure 7, which shows the damping rate in one 
period of each wave mode. Figure 7 (top panel) shows that 
the Alfven-whistler modes are abruptly damped at kp, ~ 1 . At 
smaller scales the most oblique modes are the least damped. 
The linear damping rate may thus play a role of a “filter that 
lets pass” only very oblique modes at small scales. This may 
explain the high oblique modes frequently observed in the SW 
(Sahraoui et al. 2010b; Naritaet al. 201 1; F. Sahraoui et al. 2011, 
in preparation). 

A similar extension of the fast mode solutions to frequencies 
higher than a> c i is shown in Figure 8. We observe frequency gaps 
in the dispersion curves caused by the ion Bernstein modes, 
known to develop in hot plasmas at quasi-perpendicular angles 
of propagation (e.g., Li & Flabbal 2001). The more oblique 
Bernstein modes are the least damped as shown in Figure 7 
(lower panel). This result clearly invalidates the magnetosonic 
dispersion u> ~ kC s of the classical-whistler mode at co > co c \ 
found in the previous section from the hot two-fluid theory. 

Figure 7 shows clearly that, while the Bernstein and the 
fast modes are less damped than the Alfven-whistler modes 
near kp t ~ 1, both the Bernstein and the fast magnetosonic 
modes become heavily damped at kpt > 3. The Alfven-whistler 


4 This solution has been used in Sahraoui et al. (2009) to interpret SW 
observations. 


1.0 


0.5 

LJ 5 ' 0.0 
cn 
< 

-0.5 
-1.0 

0.01 0.10 1.00 10.00 0.01 0.10 1.00 10.00 
kp. kp. 

Figure 4. Phase of the electric field component E y of the Alfven-whistler and the fast modes for the angles of propagation 80° sf Urb ^ 89° : it is positive (negative) 
for right (left) hand polarized waves. 
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Figure 6. Dispersion relations of the Alfven- whistler modes (blue) and their damping rates (red) for the angles of propagation 80° ^ #kB ^ 89? 99. 

(A color version of this figure is available in the online journal.) 



Figure 7. Exponential of the damping rates of the Alfven-whistler, fast, and Bernstein modes, normalized to the frequency of each wave for angles of propagation 
80° < 6> kB < 89°. 


modes appear to be the less damped ones at scales kp, kp 3. 
Furthermore, among these modes the most oblique one is the 
least damped, which is the KAW mode at 0 kB = 89° 99 in 
Figure 7. 

In Figure 9 we tested the validity of these results for different 
high Pi values and with 7’ = T e . We see clearly that the 


conclusions above remain valid and that the ratio 7) / T e does 
not modify much the results except for slight changes of the 
dispersion curves near the harmonics of ions. This can be seen 
for instance by comparing Figure 9 (middle panel) to Figure 3 
which was plotted for the same Pi = 2.9 but with 7) = 57). 
It is important, however, to note that when p t is becoming 
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Figure 8. Dispersion relations of the fast and the ion Bernstein modes (blue) and their damping rates (red) for the angles of propagation 80° ^ #kB ^ 89°. 

(A color version of this figure is available in the online journal.) 


smaller (i.e., decreasing toward 1) the damping of Alfven- 
whistler modes by cyclotron resonances become important for 
some angles of propagation, as can be seen in Figure 9 (top 
panel). We tested other values of /!,■ (not shown here) and found 
that the conclusions above remain valid. Moreover, we found 
that higher values of yS,- allow extending the Alfven-whistler 
modes to even smaller scales (for /?,- = p e = 25, the highest 
tested value, the modes were found to extend up to kp, ~ 30). 
On the contrary, for /l, ~ 1 we found that none of the Alfven- 
whistler modes (when 80° ^ 0 < 90°) can propagate at a> > u > c i 
because they are strongly damped at &> ~ <w c i. 

Now one can ask the question: which of these plasma modes, 
fast, Bernstein, KAW, or whistler, is likely to carry the energy 
cascade of turbulence down to the dissipation scales in the 
limited range of SW parameters studied here? While all these 
modes might contribute to the energy cascade in the SW, based 
on linear damping rates and assuming that (quasi-)linear theory 
is applicable to small-scale SW fluctuations, we can conclude 
that the KAW branch ( u> < tu c i) is more likely to be observed 
in the data than the whistler branch (a> > pj ci ). However, 
damping rates alone cannot rule out the presence of fast or 
ion Bernstein modes in the transition range 0.5 ^ kp, ^ 3 
where they are found to be less damped than the Alfven-whistler 
modes. Why, then, rule out the fast and the Bernstein modes in 
this range of scales? It has been often argued (e.g., Howes 
2009) that the Bernstein modes are highly electrostatic and 
therefore cannot account for the magnetic spectra observed in 
the dispersion range in the SW. This is not totally correct as 
can be seen in Figure 10, which shows the three components 


of the electric field of each mode. With k chosen in the XZ 
plane, and considering the high oblique angles studied here 
(80° ^ 0 kB 89 ), the Ex component (in red in the plot) 
represents essentially the electrostatic field. It is instructive to 
observe that the electromagnetic part of the Bernstein mode 
dominates over the electrostatic part at scales 0.2 ^ kp, ^ 1 
and that the electrostatic component dominates only at smaller 
scales. More interestingly, we observe that the Alfven-whistler 
modes are essentially electrostatic at all scales. This, in fact, 
is not surprising: in the MHD limit the field perturbations of 
the shear Alfven mode <5E and 5B and the mean field Bo are 
orthogonal to each other. At quasi-perpendicular angles the 
electric field <5E is therefore necessarily quasi-longitudinal (i.e., 
parallel to k). 

Figure 10 shows that only the fast mode has a dominant 
electromagnetic field up to scales kp ,■ > 3 where it becomes 
comparable to the electrostatic part. The three modes have very 
small parallel electric field (given by the E z component in 
the plot). Therefore, Figure 10 shows clearly that one cannot 
rule out any mode only on the basis of the strength of its 
electrostatic component. Instead we suggest using the magnetic 
compressibility (Song et al. 1994; Lacombe et al. 1995), in 
addition to the damping rates discussed above, to rule out the 
fast and the Bernstein modes in SW observations. 

We define the magnetic compressibility by Cg = SBf/SB 2 
(where SB 2 is the total magnetic power). It is a well-established 
result that SW magnetic compressibility is small in the inertial 
range, i.e., Cg ~ 0.1, then it increases to Cg ~ 0.4 in the tran- 
sition and the dispersive ranges (Sahraoui et al. 2010b; Kiyani 
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kp, 


Figure 9. Dispersion relations of the Alfven-whistler modes (blue) and their damping rates (red) for the angles of propagation 80° ^ #kB ^ 89° and for different 
values of (all have 7} = T e ). 

(A color version of this figure is available in the online journal.) 


et al. 2011). In Figure 1 1 we plot the magnetic compressibilities 
of the three modes for the same angles of propagation as be- 
fore. The magnetic compressibility profile of each mode shows 
clearly that the Alfven-whistler modes fit SW observations bet- 
ter, at least at the scales kp, > 1, since both the Bernstein and 
the fast modes show the dominance of parallel power over per- 
pendicular power ( Cb ~ 1). However, it is important to note in 
Figure 1 1 (top panel) that all of the Alfven-whistler modes (at 
different angles) have the same profile of magnetic compress- 
ibility. This implies that, at the oblique angles studied here, 
one cannot answer the question as to which branch, KAW or 
the whistler, dominates in the data solely from the measure- 
ment of the magnetic compressibility. This contrasts with the 


conclusions given in Gary & Smith (2009). As discussed above, 
the computed damping rates showed that the whistler branch is 
more damped than the KAW branch at small scales. 

4. DISCUSSION 

The previous results show that in high /3, plasma and at 
high oblique angles of propagation, there is one solution that 
extends from the MHD to the electron scales: it is the Alfven- 
whistler mode. The whistler mode discussed here develops as a 
continuation at high frequency of the shear Alfven mode and not 
of the fast magnetosonic mode. This result rises some exceptions 
to the general claim made by Howes et al. (2011a) in the absence 
of whistler modes in small-scale SW turbulence because of the 
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Figure 10. Electric field components E x (red), E y (green), and E z (blue) of the Alfven-whistler, fast magnetosonic, and the ion Bernstein modes normalized to their 
total electric field, computed for the angles of propagation 80° ^ ^ 89° . The E x component (red) is essentially the electrostatic part of the electric field in each 

mode, while E z (blue) is the parallel component (to Bo) of the electric field. 

(A color version of this figure is available in the online journal.) 




Figure 11. Magnetic compressibility of the Alfven-whistler, fast, and Bernstein modes for the angles of propagation 80° 89°. 


absence of fast magnetosonic modes in the inertial range. The 
fast mode strongly resonates with ions near o> cl and splits into 
different ion Bernstein modes at higher frequencies. In contrast, 
the shear Alfven mode extends at scales kp ,■ > 1 to frequencies 


either larger or smaller than Wa, depending on the ratio k\\/k±. 
The same mode can thus be called whistler (o> > oj c \) or KAW 
(o> < a> c i) depending on this ratio. If k\\/k± < then the 
Alfven mode follows the KAW mode with a> < a > c i at all scales 
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(even at kp e > 1). For larger values of k\\/k±, the Alfven mode 
develops a whistler branch with io > to cl at a given scale, but 
the damping rate increases as well. Based on the damping rates 
one would expect to observe more oblique KAW than oblique 
whistlers. This result may be considered as another, rather 
simpler, alternative to explain the strong anisotropies (k± 'p> k0 
observed in the SW. Indeed, there are several predictions from 
different MHD and kinetic turbulence models, e.g., critical 
balance conjecture (Goldreich & Sridhar 1995), that predict 
stronger anisotropy at small scales (Shebalin et al. 1983). Here 
we argue that the linear kinetic damping of the Alfven- whistler 
modes as discussed above may as well explain such observed 
strong anisotropy: the more oblique Alfven-whistler modes are 
indeed the less damped ones. 

From this result it appears clear that the question as to which 
branch, KAW or whistler, dominates in the data reduces to the 
question as to how oblique is the wave vector k. The answer 
requires accurate measurement of the angle which in turn 
requires simultaneous measurement of k\\ and k±. This is now 
possible using the A-fi lteri ng technique on Cluster spacecraft 
data (Sahraoui et al. 2006, 2010b; Narita et al. 2010a). However, 
as one can see from Figure 3, a change in the angle of 
propagation from Out = 80° to 0 k n ~ 90° will change the 
physics from the whistler branch to the KAW branch. Measuring 
the angle 0kB with an accuracy better than 10° is quite difficult, 
even with the powerful multi-spacecraft techniques, owing to the 
various sources of uncertainties in the measurements (Sahraoui 
et al. 2003b, 2010c). 

When simultaneous measurement of k\\ and k± is not possible, 
for instance when only single spacecraft data are available, 
another alternative, evoked briefly in Sahraoui et al. (2009), 
exists to test which mode (whistler or KAW) is present in the 
data. As we will show below, this method is however limited to 
situations where turbulence is not confined to the very oblique 
angles studied here. Let us discuss it in more detail here. The 
general formula relating the frequency of the wave in the SW 
rest frame to to the measured one onboard the spacecraft <w sc is 
given by 

co sc — u> + k Vs W cos 0kv> (14) 

where V sw is the SW speed forming an angle (\v with the k 
vector. If the KAW branch ( to < a> c \) dominates in the dispersive 
range down to the electron scale, then the Taylor frozen-in- 
flow assumption should be valid, meaning that Equation (14) 
reduces to a> sc ~ k 14, w cos 0 k \. Now, let us assume furthermore 
that two breakpoints occur in the turbulence energy spectra at 
the ion and the electron gyroscales pi and p e (Sahraoui et al. 
2009). As shown in Figure 12, these breakpoints should be 
observed in the spacecraft frame, respectively, at the Doppler- 
shifted frequencies f Pi = tOpJln = F sw cos0kv/27rp; and 
fp € = co p J2n = V sw co&O^y /2n p e . Assuming a similar angle 
0kv at the ion and the electron scale, the ratio between the two 
frequencies should fulfill the relation 



From Figure 1 2 one obtains the ratio f Pe / f Pi ~ 42 considering 
that 7j = T e . In the data studied in Sahraoui et al. (2009) 
Equation (15) has been accurately verified: with T,/T e ~ 5 the 
predicted ratio from Equation (15) is f Pe /f Pi ~ 95, which was 
very close to the ratio between the frequencies corresponding 
to the observed electron breakpoint at f Pe ~ 40 Hz and ion 
breakpoint at f Pi ~ 0.4 Hz. 



Figure 12. Frequency signatures of the ion and electron gyroscales a frozen-in 
flow approximation (red) and two whistler modes at two different angles of 
propagation (7) = T e ). The indicated frequencies are given in the text. 

(A color version of this figure is available in the online journal.) 

If, however, relation (15) is not fulfilled by two observed 
breakpoints, then this would imply that the Taylor assumption 
fails, or equivalently, that a dispersive branch of the whistler 
type exists within the data. To show that let us now assume 
that the whistler branch is dominant in the dispersive range. In 
this case the Taylor hypothesis will naturally fail (since phase 
speeds comparable to V sw are present). At the very small scale, 
i.e., near p e , one can even expect that the phase speed of the 
whistler mode, being linear in k (since to oc k 2 ), will be larger 
than F sw as can be seen in Figure 12. In this case Equation (14) 
yields &> sc ~ to. This means that the relevant temporal scale will 
be given, not by “Doppler-shifting” p e , but by the dispersion 
curve of the whistler mode in the plasma rest frame. From 
Figure 12 one can see that the scale p e (—d e ) will yield a 
different breakpoint frequency depending on the angle 0 k \\'. 
to i ~ 500&) c i and — 1000&> c i, respectively, for f?kB = 60° and 
$kB = 10°. With fa ~ 0.5 Hz, the breakpoint at the electron 
scale should thus be observed onboard the spacecraft at the 
frequencies 500 Hz and 250 Hz, respectively. These frequencies 
are much higher than those reported in Sahraoui et al. (2009), 
which validates the conclusions of that work that turbulence 
was more consistent with the KAW modes (or quasi-stationary 
fluctuations) than with whistlers modes. At higher oblique 
angles (0 k u > 60°) one can extrapolate from Figure 12 that the 
phase speed of the whistler branch will become comparable or 
smaller than V sw , therefore the Doppler shift cannot be neglected 
in Equation (14). An accurate estimation of both terms in the 
right-hand side of Equation ( 14) is then necessary. As mentioned 
above, this shows the limit of this simple method and emphasizes 
again the difficulty in distinguishing between the two branches, 
KAW and whistler, of the same mode at the very high oblique 
angles. The estimations given here can be tested on S W magnetic 
power spectra that can be measured by the high-time resolution 
Cluster search coils (up to 225 Hz in the spacecraft frame of 
reference Sahraoui et al. 2009, 2010c) to search for data intervals 
when parallel or moderate oblique whistler mode turbulence 
may exist (i.e., intervals not fulfilling Equation (15)). However, 
very often, the level of the magnetic turbulence in the SW is 
not sufficiently high and can hit the sensitivity floor of the 
magnetometers at frequencies as low as 20 Hz. This limitation 
emphasizes the need in the future space missions for new search 
coils that have higher levels of sensitivity in the dispersive and 
the dissipation ranges of SW turbulence. 
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5. CONCLUSIONS 

Using the hot two-fluid and the Vlasov theories we studied 
key properties of the linear plasma modes under realistic SW 
conditions at 1 AU, namely, Pi > A ~ 1, Ti > T e and focused 
on high oblique angles of propagation 89° ^ (9kB < 90° as 
observed frequently from the Cluster multi-spacecraft data. We 
discussed the relevance of each of the main modes (the KAW, 
the whistler, the fast magnetosonic, and the ion Bernstein) to 
carry the energy cascade below the ion gyroscale (the slow 
magnetosonic mode being heavily damped by kinetic effects). 
We showed in particular that the whistler branch develops as a 
continuation at high frequency of the classical (shear) Alfven 
wave known in the MHD limit, and not of the fast magnetosonic 
mode. The Alfven mode extends indeed at scales kpi > 1 to 
frequencies either larger or smaller than o> cl , depending on 
the anisotropy k^/k±. The same mode is thus called whistler 
(<y > <w c i) or KAW (a> < co c j) depending on the anisotropy ratio. 

Unlike the two-fluid model, the Vlasov theory shows that the 
fast magnetosonic mode undergoes strong cyclotron damping 
at co ~ &> c i and splits up into ion Bernstein modes at to > a > CI , 
in agreement with the analysis of the wave polarization. We 
also showed that ruling out the fast and the Bernstein modes 
from the cascade process in the dispersive range based on their 
electrostatic nature is incorrect, since the Alfven- whistler modes 
are shown to be electrostatic as well at high oblique angles. We 
suggest rather that one must combine magnetic compressibility 
and damping rates to rule out those modes: both were shown to 
be highly compressible (i.e., B 1 Bj) in disagreement with 
typical SW observations of magnetic field data. In addition, 
those modes were shown to be heavily damped in the dispersive 
range as compared with the Alfven-whistler mode. We showed, 
finally, that the modes on the whistler branch are more damped 
than those of the KAW branch, which makes the latter a more 
relevant candidate to carry the energy cascade down to electron 
scales. We also demonstrate a way to test the presence of either 
branch in SW spacecraft data and discussed the limitation of the 
existing spacecraft observations that need to be kept in mind in 
the planning of the future space missions. 

Finally, it is important to keep in mind that the conclusions of 
this work may not be applicable to all SW observations because 
the present study is fully linear (one can question the relevance 
of linear theories to SW turbulence, see, for instance, Dmitruk 
& Matthaeus 2009; Howes 2009) and restricted to the range 
of parameters discussed above. Moreover, the role that can be 
played by several plasma instabilities (both at ion and electron 
scales) as reported in several space observations (Sahraoui et al. 
2006; Hellinger et al. 2006; Bale et al. 2009) is not addressed 
here. 

This work is part of the project THESOW funded by F Agence 
Nationale de la Recherche (ANR, France). F. S. thanks G. Howes 
for stimulating discussions. 
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